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In a recent paper, we combined the technique of bosonization with the concept of a Rayleigh 
dissipation function to develop a model for resistances in one-dimensional systems of interacting 
spinless electrons arXiv:1011.5058 . We also studied the conductance of a system of three wires 
by using a current sphtting matrix M at the junction. In this paper we extend our earlier work 
in several ways. The power dissipated in a three-wire system is calculated as a function of M and 
the voltages applied in the leads. By combining two junctions of three wires, we examine a system 
consisting of two parallel resistances. We study the conductance of this system as a function of the 
M matrices and the two resistances; we find that the total resistance is generally quite different 
from what one expects for a classical system of parallel resistances. We will do a sum over paths 
to compute the conductance of this system when one of the two resistances is taken to be infinitely 
large. Finally we study the conductance of a three-wire system of interacting spin-1/2 electrons, 
and show that the charge and spin conductances can generally be different from each other. 

PACS numbers: 73.23.-b, 73.63.Nm, 71.10.Pm 



I. INTRODUCTION 



For non-interacting electrons, the conductance of a bal- 
listic quantum wire is well-known to be quantized in units 
of 2e^ /h at low temperatures This remains valid 

when interactions between the electrons in the wire arc 
taken into account, provided that there are no sources 
of backscattering (such as impurities and junctions) and 
that the wire is connected to leads where there are no 
interactions 0-111 • Thus, if the wire is modeled as a 
Tomonaga-Luttinger liquid (TLL) and the interaction 
strength is given by the Luttinger parameter Kwi the 
conductance of a clean wire does not depend on Kw- 
This breaks down if there are impurities in a wire with 
interacting electrons. In that case, the impurity strengths 
vary with the length scale according to some renormaliza- 
tion group (RG) equations, and the conductance depends 
on Kw and other parameters like the length of the wire, 
the distances between the impurities, and the tempera- 
ture considerable amount of work has also been 
done on junctions of several quantum wires both theo- 
retically [15| - |23 | and experimentally [lH, [2^. In these 
systems, the conductance matrix again becomes length 
scale dependent due to the interactions. A junction of 
three quantum wires with interacting spin-1/2 electrons 
has been studied in Ref. [l^] , and it has been found that 
some of the fixed points of the RG equations exhibit dif- 
ferent charge and spin conductances. Thus impurities 
in a single wire or a junction of three or more wires ef- 
fectively give rise to a resistance which leads to power 
dissipation. There have been some studies of power dis- 
sipation on the edges of a quantum Hall system and 
at a junction of quantum wires [T^ . However, there has 
been relatively little discussion in the literature of the 
effects of an extended region of dissipation (a resistive 
patch) within the framework of TLL theory or bosoniza- 
tion which is the most efficient way to study the effects 
of interactions [isll. 



To remedy this situation, we recently introduced a for- 
malism which can combine the technique of bosonization 
with the classical notion of resistance; both a single wire 
and a system of three wires with a junction were studied 
using this formalism j^. Our analysis was restricted to 
spinless electrons and zero temperature. In contrast to 
this, Ref. dl] considered the effect of an extended region 
of inhomogeneity in a quantum wire at low temperatures; 
it was shown that this leads to weak backscattering which 
gives rise to a resistance which is linear in the tempera- 
ture. 

It is useful to briefly recapitulate our earlier work [ij] . 
We introduced the resistances phenomenologically using 
a Rayleigh dissipative function [2^. Our treatment was 
classical in the sense that the resistance was taken to be 
purely a source of power dissipation, and the microscopic 
quantum mechanical origins of the resistance were not 
specified. This is equivalent to treating scattering by the 
resistance as a phase incoherent process; a consequence 
of this is that the resistances of different patches add up 
in series with no effects of interference. We found expres- 
sions for the conductance of a single wire and of a three- 
wire system with a junction (this was described by a cur- 
rent splitting matrix M which is orthogonal) in terms 
of the resistances in the wires. The conductance can be 
calculated even when the Luttinger parameter Kw, the 
velocity vw of the quasiparticles and the resistivity r all 
vary with the spatial coordinate x in an arbitrary way 
in the different wires. Remarkably, we found that the 
conductance of the three-wire system is independent of 
the Luttinger parameter Kw (which is determined by the 
strength of the interaction between the electrons near the 
junction) if the matrix M is invariant under time rever- 
sal; this remains true no matter what the values of the 
resistances are in the three wires. 

In this paper, we will extend the ideas introduced in 
Ref. [iJl in several ways. For the sake of completeness, 
we will briefly present our earlier results for a single wire 



in Sec. II and a three-wire system in Sec. III. We will 
then examine the issue of power dissipation in a three- 
wire system in Sec. IV. The dissipation occurs due to 
the resistances in the wires and the contact resistance in 
the leads. We will study the dependence of the power 
dissipation on the junction matrix M and the relative 
magnitudes of the voltages applied in the leads of the 
system. In Sec. V, we will study a system of two resis- 
tances in parallel; the system consists of two junctions 
of three wires. We will make a detailed study of the 
total conductance in terms of the two junction matri- 
ces M and the two resistances; the conductance will be 
calculated using a scattering approach. We find that the 
classical expression for the effective resistance of two par- 
allel resistances is recovered only in one special case. In 
Sec. VI, we will carry out a sum over paths to calculate 
the conductance of the same system in the special case 
that one of the two resistances is infinitely large. In Sec. 
VII, we will generalize our results to the case of spin-1/2 
electrons. We will argue there that the charge and spin 
conductances can generally be different from each other. 
We will summarize our results in Sec. VIII. 



II. SINGLE WIRE 



which we will take to be — L/2 < a; < L/2. (For non- 
interacting electrons, K = \^ while for short-range re- 
pulsive interactions like a screened Coulomb interaction, 
we have K < 1). The limit of noninteracting electrons, 
with K = 1 and v = vp (the Fermi velocity), is used to 
model the two- or three-dimensional Fermi liquid leads 
situated in the regions |a;| > L/2. In the leads, the fre- 
quency and wave number of a plane wave are related by 
uj = VF\k\. The electron charge density n and current j 
can be expressed in terms of the bosonic field by the rela- 
tions: n = —edx<p/V^ E^nd j — edt4>/y/Tr, where e is the 
electron charge; these satisfy the equation of continuity 
dtn + dxj = 0. To describe resistances phenomenologi- 
cally, we introduce a Rayleigh dissipation function 



J" 



dx r j , 



(2) 



where the resistivity r can vary with x. We then obtain 
the following equation of motion as described in Ref. [l^l 
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Note that we have set ft = 1, so that e^/(27r) = /h. 




FIG. 1. (i) Schematic picture of a resistive region connected to 
Fermi liquid leads (FLL). (ii) A grapii showing the resistivity 
as a function of the length; dx r{x) — R. 



A. Equation of motion 

We start with the analysis of a single wire containing 
interacting spinless electrons. In the absence of any im- 
purities, the bosonic Lagrangian for the system is given 

by 



dx 
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where K and v respectively denote the Luttinger param- 
eter and velocity of the bosonic quasiparticles; we will 
allow K and v to vary with x within a finite region 



B. Scattering approach 

We will now derive an expression for the DC conduc- 
tance for a general resistance profile r{x), and will show 
that only the total resistance R = J dxr{x) of all the 
resistive patches appears in the final expression for Udc- 
Also, we will show that a 5-function resistance with the 
same integrated value of R gives the same value of the 
conductance. 

As mentioned earlier, we set K = 1 and v = vp i'n the 
leads. However, K{x) and v{x) can have any profile in the 
wire region given by |a;| < a. Similarly, we will assume 
that the resistivity r{x) = for \x\ > a, but can have any 
profile in the region |x| < a such that /"^ dx r{x) = i? 
[one such profile is illustrated in Fig.[T](ii)]. As described 
in Ref. [ij], in the scattering approach, a plane wave 
with frequency uj is incident on the resistive patch, the 
reflection and transmission amplitudes arc calculated as 
functions of w, and finally the limit w — > 0^ is taken to 
obtain the expression for adc- For a plane wave incident 
from left with k = lu/vf, the spatial part of the solution 
(j)k{x,t) = /fc(a;)e~*"* outside the resistive patch is given 
by 



/fe = e*^" + Sfc e-*'^'" for x < -a, 
= tk e^^^ for a < X. 



(4) 



Up to zero-th order in ui (and hence in k) the solution of 
Eq. (|3]) is = /fc = c, where c is a constant. Therefore 
we can see from Eq. ^ that 



Sk 



up to zero — th order in u. (5) 
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Since we are eventually interested in the limit w — >■ 0"*", 
we can work out the solution to Eq. ([3|) up to first order 
in uj and then take the limit 0+. We first rewrite 

Eq. ([3]) up to first order in w, 

2 

- {^d.h) - r h ^ Q. (6) 

In the second term in Eq. ([5]), we can replace fk by a 
constant c (= tk up to zero-th order in to) since this term 
has a factor of uj already. Integrating this equation from 
X = —a — e to a + e, we obtain ikvpitk ~ (1 ^ Sfc)] = 
—iujtk{e'^R/TT) which is the same as 

+ + = 1. (7) 

Solving Eqs. ^ and we obtain = 1/[1 + 
Note that this expression for tk is correct only in the 
limit a; — > 0+ (or A: — >■ 0+). Hence, 

e2 

ZTT 



One can easily redo the calculations from Eq. Q to 
Eq. ([5]) and see that the result remains unchanged if we 
choose a (5-function resistivity profile described by r(x) — 
RS{x — Xq), where |xo| < a. We will use this fact later 
in Sec. I VII We thus see that adc does not depend on 
the precise functional forms of K{x), v{x) and r{x), as 
long as K and v are equal to the constants 1 and vp in 
the leads, and R is the total resistance of the wire. An 
implication of this is that the effective resistance of two 
or more resistive patches in series will be given by the 
sum of the individual resistances. 

In this context, it is worthwhile to look at the phase 
coherent and phase incoherent transport in the literature 
(see for instance pp. 125-129 of Ref. (H). In general, 
the effective resistance of two resistances can have an 
extra term that depends on some phase factors at the two 
resistive patches in addition to the sum of the individual 
resistances. However, the effective resistance reduces to 
the sum of individual resistances in the incoherent limit. 
In this sense, our formalism assumes that the system is 
phase incoherent. 

III. THREE WIRES WITH A JUNCTION 

Let us consider a junction of three TLL wires as shown 
in Fig. [5] We will assume that each wire has three re- 
gions: 

(i) < < Lii 7^ — the wire region around the 
junction where K{xi) = Kw and w(.t,;) = vw', elsewhere 
K{xi) = 1 and v{xi) — vp, 

(ii) Lii < Xi < Li2 — the dissipative region where 
r{xi) = Tio; elsewhere r{xi) = 0, and 



(iii) Xi > Li2 — the semi-infinite leads. 

Here i labels the wires, and on wire i, the coordinate 
Xi runs from to oo, with Xi ~ denoting the junction 
point. The regions Xi > Li2 model the two- or three- 
dimensional leads which are assumed to be Fermi liquids 
with no interactions between the electrons; we therefore 
set A' = 1 and v ~ vp in those regions. 




FIG. 2. Schematic diagram of a three-wire junction with in- 
teracting regions close to the junction (pink solid region), dis- 
sipative regions (green zig-zag region further from the junc- 
tion), and Fermi liquid leads (brown shaded region furthest 
from the junction). 



A. Green's function approach 

We now follow Ref. Q and write 

J ^e-*'^V,j,„(a;,, x;)£;^(x;), 

(9) 

in the linear response regime, where Ei^{x'j) is the Fourier 
component of the electric field E{x'j,t) on wire j, and 
o'ij.ui{xi, x'j) is the nonlocal ac conductance matrix. Using 
the Kubo formula, we then obtain 

9 - 

where uj = —ico, and 

Q,,,^{x,,x'^) = —{T;U^,,T)<p,{x'^,Q))e'^^ (11) 

is the Fourier transform of the bosonic field in imaginary 
time, T — it. The expressions for t and w in terms of t 
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and uj arise as follows. The actions in real (Minkowski) 
time t and imaginary (Euclidean) time r are given by 

Sm = j J dtdx[^{Mf - ^{d.<P)\ 

Se^J J drdx [^{dr^f + 2^(5-<^)'] (12) 

respectively. The requirement that the exponentials ap- 
pearing in a path integral formulation be equal to each 
other, i.e., e*'^*^ = e^'^^, implies that r = it. Secondly, 
we want an outgoing plane wave in the lead of wire i to 
be given by e'("^»/"^~"*) in terms of a real frequency w 
and Q-'^^i/^F-^^'^ in terms of an imaginary frequency w 
(this expression for an outgoing wave will be used in the 
boundary condition (iii) given below). The two planes 
will be identical if w = —iu. 

The Green's function satisfies the equation 

/ v{xi) \ Gp- e^Cd 

Gij,cjixi,x'j) = Sij S{xi - x'j), (13) 

with the following three boundary conditions: 

(i) Qij.u{xi, x'j) is continuous at Xi — x'^ (where < < 

Lj2) and -;^^9a..^jj,o(a;i,a;j-)I^Ce = 

(ii) Gij,ui{xi,x[j) and -^^^^dx,Gij,a{x.i,x'^) are continu- 
ous at Xi = Lii and Li2, 

(iii) if g^j,a{x,,x'j) = Ay-e'^^-/""' + Bye-""-/""' for < 
Xi < min(a;^, Lii) (5^ -I- Lii{\ — 5ij), then B = —M A, 
where M is the current splitting matrix at the junc- 
tion [H-lil. 

The boundary condition in (iii) encodes the fact that the 
incoming and outgoing currents (and hence the bosonic 
fields) at the junction are related by the matrix M . Var- 
ious constraints at the junction such as current conser- 
vation and unitarity of the evolution of the system in 
real time (i.e., no power is dissipated exactly at the junc- 
tion) imply that each row and column of M must add 
up to unity and that M must be orthogonal. It turns 
out that for a junction of three wires, the possible M 
matrices must belong to one of two classes both of which 
are parameterized by a single parameter 9 (l5l - [2^ : (a) 
det(Mi) = 1 and (b) det(M2) = —1; these can be ex- 
pressed as: 
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(14) 



where a = (1 + 2cos6')/3, b = (1 - cos0 + A/3sin6')/3, 
and c = (l-cos6'-^/3sin6')/3. We note that {M2Y = 1 
for any value of 9] this relation will be used below. 

Note that by introducing the orthogonal matrix M, we 
have made the assumption that there is no dissipation 
exactly at the junction. The analysis is simpler if the 
junction, which governs how the incoming currents are 
distributed amongst the different wires, is separated from 



the dissipative regions which lie some distance away from 
the junction. 

Solving Eq. (jl3p with the above boundary conditions 
and taking the limit w — > 0"^, we obtain the following 
expression for the dc conductance matrix: 

G=- [1 H- M + Kw{t - M){t + -R)]-i 

X [1-A/], (15) 

where IR is a 3 x 3 diagonal matrix with = Ri = 
i'io{Li2 — Lii); we note that Ri is simply the total re- 
sistance in wire i. The conductance matrix relates the 
outgoing current Ii to the potential Vi applied in lead i 
as Ii = J2j C!ijVj. One can show in general that each 
row and column of G must add up to zero; the columns 
adding up to zero is a consequence of current conserva- 
tion (^j Ii must be zero), while the rows must add up to 
zero because each of the Ii must vanish if the V^ 's have 
the same values in all the wires. 

Eq. p5|) can be understood as a combination of the 
conductance of a system with no resistances {Ri — 0) 
and resistances Ri on the three wires. Let us denote the 
conductance with no resistances by 

Go = - [1 + M + Kwit-M)]-^ [1-AI]. (16) 

TT 

If V/ denote the potentials at the points Xi ~ Ln (i.e., 
the points which lie after the interacting regions but be- 
fore the dissipative regions), we have Ii = J2jGo,ijVj- 
Further, /; = {V^ — Vi)/Ri. Combining these equations 
with /; = GijVj, we obtain 

G = [1 - GoR]-' Go. (17) 
This relation will be used in Sec. IVIII below. 

B. Scattering approach 

Eq. (|15[) can be derived in general using the equation of 
motion approach in the limit a; — > 0"*" in the same way as 
described above for the single wire case. Let (j)i{xi,t) = 
/i(xi)e-*"* denote the bosonic field in real time which 
satisfies Eq. ([3]) in wire i. Current conservation at the 
junction implies that '^^dt4>i{xi — > -t- e, i) = which 
implies that J^i fi{xi — > + e) = 0. [Note that this is a 
different condition than the one used in the single wire 
case where f{x) was assumed to be continuous every- 
where; for the case of more than two wires, it is more 
convenient to assume J^i fi {xi ^ + e) = rather than 
the continuity of fi {xi ^ + e) between different values 
of i] . We now assume that 

/.(x, ^0 + e) = a,e-^'='"- + fe.e''^'"-, 

/,(x, ->oo) = a,e-''^"- + Ae^'^"', (18) 

where a^, ai denote the incoming fields and 6^, f3i de- 
note the outgoing fields. Assuming that the Luttinger 
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parameter and the velocity are given by Kw and vyy as 
— > + e and by 1 and vp as Xi oo, we must have 
vwk' = vpk ~ Ld. In the limit lu ^ 0^, the field on 
wire i is given by + 6^ at Xi — >■ + e and by ai + /3i 
as Xi — >■ oo at zero-th order in w, k and k'. Since fi{xi) 
equal to a constant is a solution of Eq. ([3|) for cj = 0, 
we must have ai + hi = ai + Next, the coefficients 
ai and bi arc related by the current splitting matrix at 
the junction, b = —M a. The coefficients ai and /3i 
must be related by the conductance matrix in the leads, 
/? = — [(27r/e^)G + D] a; this follows from the statement 
that J2ii^i+ Pi) = + = 0- Finally, we integrate 

Eq. ([3]) from + e to oo, ignoring the first term 

which is of order and setting the third term equal 
to ~iuj{e'^/TT){ai + Pi)Ru, where Ri^ = Jj^'^ dxtr{x,). 
This gives the equation j3i — ai — {1/Kw){bi — ai) = 
— (e^/7r)(ai + /3i)IRii, where we have used the fact that 
v-^k' = Vpk = u and taken the limit w — >■ 0^. Using all 
these equations, we recover Eq. ([TS]). We thus see that 
the precise profiles of K{xi), v{xi) and r{xi) in the differ- 
ent wires are not important; all that matters is that the 
values of K and v are given by Kw and vw as cc^ — > + e 
and by 1 and vp as Xi — >■ oo, and that Ra — J dxir{xi). 



C. Conductance for the Afi and AI2 classes 

Within the Mi class, the case = is trivial because 
Mi(0) = 1 and G = 0. We now consider all other values 
of 6. We find that in general G depends on K^y, 9, and 
the resistances Ri = rio{Li2 — Ln). [An exception arises 
for the case = n where we find that G is independent of 
and depends only on the Ri. As shown below, this 
occurs whenever M"^ = 1 which is true for Mi(7r) and 
also for the M2 class for any 9.] The dependence of G 
on Kw for the Mi class is to be contrasted to the case 
of a single wire where the conductance is independent of 

Kw m- 

In the AI2 class, we find that although G depends on 
6 and the Ri, it is completely independent of Kw- In 
Eq. (Hll), we write 1 + M + Kw{t - M)[l + ie'^/Tr)R\ = 
A + B, where A = 1 + M + Kw{t - M) and B = 
{e"^ /TT)Kw{t — M)R. We can then use the relations that 
AP = 1, = [Kw{^ + M) + (1 - M)]/{'iKw) and 

{A + = - A-^BA-^ + A-^BA-^BA-^ 

to show that G does not depend on Kw for any choice 
of 9 and Ri. The exact expression for G turns out to be 

3(1 - M2) 
^-'^ D ' 
where D = 2{gi + g2 + ga) + cos 9{gi + g2 - 2^)3) 

^V3sm9{gi - g2), (19) 

where = 1 + (e^/7r)i?,;. We thus see that G docs not 
depend on Kw, the Luttinger parameter in the wire re- 
gions. 



IV. POWER DISSIPATION 

In our model, there is no power dissipation exactly at 
the junction since the current splitting matrices Mi and 
M2 are orthogonal. Power dissipation occurs only at the 
resistive patches and in the leads due to the contact re- 
sistance. The power dissipation at the contact resistance 
occurs due to the energy relaxation of the electrons in the 
leads (reservoirs) which are maintained at some particu- 
lar chemical potentials. Classically, if a current Iq passes 
through a resistance Rq, the power dissipated Pq is given 
by Pq ~ I§Ro. For a three- wire junction, we can define 
the power dissipated in two equivalent ways as follows 



P= - ^ K 



i=l 
3 



and P= E^H^^ + i^)- 



(20a) 
(20b) 



We have verified analytically that these two definitions 
give the same result. (A minus sign appears in Eq. (j20ap 
because we have defined the li to be outgoing currents). 

We know that when the voltages applied in all the 
three leads are equal to each other, there should be no 
current in any of the three wires, and hence the power 
dissipated should be zero. To incorporate this fact, we 
choose a new coordinate system for the Vi 's known as the 
Jacobi coordinates 







Vb 




Vc 





l/\/3 1/V3 l/%/3 
1/72 -1/72 
l/\/6 i/Ve -2/V6 



(21) 



In this coordinate system, the power dissipated P does 
not depend on the voltage Vq. Further, it turns out that 
Vf, and Vc can be parameterized in such a way that P de- 
pends only on one of two parameters; this parametriza- 
tion is different for junctions described by M\ and M2 as 
we will see below. Since we know from Eq. (j20b[) that the 
power dissipated at the resistance Ri is l^Ri, for simplic- 
ity we will only consider the case of i?i = (for i = 1, 2, 3) 
for further analysis. 



A. Power dissipated for the M\ class 

In this case, the A/-matrix is invariant under a cyclic 
permutation of the wires 1, 2 and 3, and the dissipated 
power P turns out to be invariant under rotations in the 
Vf, — Yc plane. If we write VJ, = Vq cos and Vh — Vq sin 0, 
the power dissipated is proportional to Vq and does not 
depend on 0. However, P depends on the parameters 
Kw and 9 since the conductance matrix G depends on 
those two parameters. The dependence of P on Kw and 
9 is shown as a contour plot in Fig. |31 The choice 9 ~Q 
decouples the three wires at the junction, making the cur- 
rents in the three wires zero; hence, the power dissipated 
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is zero for this choice. The power dissipated is maximum 
on the contour 9 = ±7r. 
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FIG. 3. Power dissipated for a Y-junction (see Fig. [2]) shown 
as a contour plot in the 6 — Kw plane for the Afi class. The 
numbers shown indicate the power dissipated along the near- 
est contour in units of /h\ we have set Vb = 1 and Ri = 0. 



B. Power dissipated for the A/2 class 

Junctions described by Af2-class are time reversal in- 
variant (the M- matrix is symmetric). In this case, the 
power dissipated is a constant over straight hnes in 
the Vb — Vc plane. If we define the variables V± = 
Vb cos{9/2) ± VcSm{9/2), then the power dissipated P 
is found to be independent of V+ and is given by P = 
(e^//i) V3- Further, P does not depend on Kw in this 
case. 



V. RESISTANCES IN PARALLEL 

At the end of SecHHwe saw that the effective resistance 
of two or more resistances in series is the sum of the indi- 
vidual resistances in agreement with the classical result. 
By the word 'classical' we mean the result obtained for 
the effective resistance by using Kirchoff's circuit laws. 
In this spirit, it is interesting to study whether the ef- 
fective resistance of two resistances in parallel will agree 
with the classical result for a similar system. 

To begin with, let us consider a general model illus- 
trated in Fig. m The three wires at each junction are 
labeled by the index i = 1,2,3, and the two junctions 
are described by current splitting matrices Mc and AI-ji. 




FIG. 4. Schematic diagram of resistances R2 and Rs in par- 
allel attached to Fermi liquid leads through junctions £ and 
7^. 



The coordinate Xi runs from 

(i) —00 to —I on the wire i = I on the left, 

(ii) — / to I on the wires i = 2, 3, 

(iii) / to 00 on the wire i = 1 on the right. 

The resistive regions on the wires i = 2, 3 lie in the range 
\xi\ < a, the total resistances being Ri. In the leads, 
K = 1 and v = vp as usual, while K = and v = vw 
in the regions a < jxij < L/2 shown in pink (solid lines). 
In the resistive regions, K[x) and v{x) can have any pro- 
files and these do not affect the final result. The incoming 
and outgoing fields at the two junctions (at xi = ±/) are 
related by the matrices M-]^ and Mc ■ A plane wave inci- 
dent from the left has the solution (j)i.k{xi, t) = fi^k e~*"* 
given by: 
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for j = 2 and 3, 
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for xi < -L/2, 

- L/2 <xi < ~l, 
(22a) 

for — I < Xi < —a, 



for 



< Xj < I 



tkC 



ikx 



(22b) 

for I <xi < L/2, 

for xi > L/2. 

(22c) 



A calculation similar to the one performed in Sec. Ill Bl 

relates (in the DC limit) the unknowns tj-p^k and Sj-p^k 
(for V = C,n and j = 2, 3) as 



(23) 



Sjv.k 

for different j's are related by the matrix M-p as follows: 





tjCk 




■ 1 


+ 9j 






tjKk 




_ Sjck _ 












SjKk 


where Qj 


= e'R,/i27r). 


The 


amplitudes tj-p^k 



w. 



where w 



w 



w. 



w. 



(out) 
V 

(out) 
C 

(in) 
C 

(out) 

n 

(in) 

n 



(in) 
V 

\-s\c.k 



Mp Wj, 



c,n, 



for V - 

hck t3C,k\ . 
[ t\c.k — S2C.k ~ S^CkV , 
[ tin,k — S2n.k — S^Ti.k^ ^ 

[-sin,k t2n,k tsTi.kf- (24) 



This is essentially same as the relation b = — M a used 
in Sec. HITB] The relations ^ and ^ along with the 
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continuity of (j)i^k and -j^dx4>i,k sX x = ±L/2 give us 
enough conditions to solve for all the unknowns, namely, 
Sk, tjV,k, Sjv.k and tfc. 

Once we know the final expression for the DC conduc- 
tance tTrfc (= ^ ifc in the limit k — >■ 0"*"), we can obtain 
the effective resistance R\ \ by subtracting out the contact 
resistance from the total resistance, 



1 



R\\ = ^ 



h 

■ 



(25) 



We have listed the results obtained for different choices 
of Mc and Mn in Table H] 



Mc{ec) 


Mn{9Tz) 


Expression for _R| 


9^0 


Ml (-61) 


R2R-j,/{R2 + Rs) 


Mi{9c) 
9c + ~9ti 




Depends on 9c-, 9ti and Kw 
as shown in Figs. [5] and [6] 


Mx(9c = 0) 


M2{9n) 


00 


M^{9c) 


M2i9n) 


Depends only on 9tz = 9 
as shown in Fig. [7] 


M2{9) 


M2{9) 


Depends only on 9 
as shown in Fig. [7] 


M2{9c + 9tz) 


M2{9n) 


00 



TABLE I. The behavior of the effective resistance for dif- 
ferent choices of Mc and M-r,. 

In an earlier paper, we observed that in the case of a 
three-wire junction, if time reversal invariancc is broken 
at the junction (i.e., if M is not a symmetric matrix), 
then the conductance matrix depends on Kw (T^l. In 
the model studied here, we find that the final DC con- 
ductance depends on only when time reversal invari- 
ance is broken at both the junctions (i.e., if both Mc and 
M-ji belong to the Mi class) and 9c 7^ —Ou- The depen- 
dence of the conductance adc on Kw {Sc/n) is shown in 
Fig. [sill]). The contour Oc ^ —d-ji = 6 corresponds 
to the maximum value of adc- On this contour, adc does 
not depend on cither Kw or 9; moreover we get an ex- 
pression for the effective resistance R\\ which agrees with 
the classical result for i?| | . 

The choice = at a junction described by Mi-matrix 
decouples all the three wires at the junction. Hence, 
the conductance of the system is zero (or equivalently 
= 00) if either of the two junctions have 9 = 0. We 
can see this in both Table U and Fig.[6l Another interest- 
ing case that results in infinite effective resistance arises 
when both Mc and Mu lie in the M2 class and 9c 9-]z- 
The choices, (i) Mc and M-ji lie in the M2 class with 
9c = (^n = and (ii) Mc lies in the Mi class and M-ji 
lies in the M2 class, with < 9c < '2-tt and 9-ji = 9, yield 
the same ^-dependent expression for the conductance, as 
shown in Fig. [7] for one particular choice of the parame- 
ters. The case Mc{9c) = M2{92) and Mc{9n) = Mi{9i) 
is redundant since we have already analyzed the case 
Mci9c) = Mi{9i) and Mci9Tz) = ^^2(6*2), and these 
two cases are equivalent. The equivalence of these two 



0.5 [- 

^ 0.4 

4-1 

° 
'a 

c 0.2 

b 

0.1 
0.0 



0.0 



0.5 



1.0 



1.5 



2.0 



FIG. 5. DC conductance is plotted versus Kw for the case 
that Mc and Mn lie in the Mi class. Parameters: 9c = 9tz = 
7r/2, £12 = 1 and £13 = 2. 











w 




==9 =^ 






K^^OTOl 


X \ \ 

0.59 0.4 


«?.2 ia.5 \ 



-3-2-10 1 2 3 



FIG. 6. DC conductance Gdc for a parallel combination of 
resistances (see Fig. [4]) shown as a contour plot in the S-r — 9c 
plane for a fixed interaction strength Kw = 0.6, when both 
Mc and Miz lie in the Mi class. We have chosen Q2 = ^, Q3 = 
2. The numbers indicate the values of adc on the contours in 
units of /h. 



cases follows from the parity transformation Xi —Xi, 
i = f,2,3, since the conductance is parity invariant for 
this system. 

Finally, as a special case, we look at a symmetric situa- 
tion treating the wires 2 and 3 on the same footing, with 
both junctions being described by the same M-matrix. 
The case M = Mi requires 9 = tt for symmetry between 
the wires 2 and 3. As 9c = —9n, for 9c = tt, we conclude 
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FIG. 7. DC conductance is plotted versus 9 for the cases (i) 
Mc and M-jz lie in the M2 class and 6c = 9n ~ 9, and (ii) Mc 
lies in the Mi class and M-r, lies in the M2 class, < 9c < 2n 
and = Parameters: ^2 = 1 and gs — 2. 



from Tabic U that i?|| = Ro^R^/iR-i + Rz)- Turning to 
the case AI = A/2, the symmetry requirement imposes 
9 = — 7r/3. The effective resistance in this case takes 
the form = {R2 + i?3)/4. An interesting implication 
of this is that when one of the resistances (say R^) is 
taken to infinity keeping the other {R2) finite, the effec- 
tive resistance R\\ also goes to infinity. [This is unlike 
the classical case where, if one of the resistances [R^) in 
a parallel geometry is taken to infinity, the effective resis- 
tance approaches the value of the other resistance (i?2)]- 
We will provide an understanding of this surprising result 
in Sec. I VII by a physically intuitive argument. 



1 Mc 2 



R2 



T 



FIG. 8. Schematic diagram of the model. The J-function 
resistivity profile is highlighted by the diamond in the middle. 
Any pulse which hits the free end of either of the wires labeled 
3 gets completely reflected. 



liquid leads as shown in Fig. [H] To illustrate the 'sum 
over paths' method, let us first consider a simpler model 
where we calculate the transmission amplitude at a In- 
junction formed by attaching a single stub to a wire as 
shown in Fig. [HI Let M-p be the M-matrix that relates 
the outgoing bosonic fields to the incoming bosonic fields 
at the junction. Of the pulse approaching the y-junction 
through the wire i, a fraction M-pji goes into the wire j 
after scattering. So the transmission amplitude for the 
path 1 2 is M-P21- Since the pulse gets completely 
reflected at the free end of wire 3, the transmission am- 
plitude for the path 1 — > 3 — > 2 is M-p23 M-p^i. 



M-p 21 



M-P23 ■ ■ Mp3i 

1 • 2 



(0 



M-P23 ■ M-psi 



(ii) 3 



FIG. 9. Transmission amplitudes for a few paths are shown 
as examples. 



VI. TRANSMISSION OF A PULSE 

Inspired by the scattering method [1, 0, [l3| , we illus- 
trate a 'sum over paths' method to calculate the DC con- 
ductance of a system. We consider a special case of the 
model studied in Sec.|V]by setting K = 1 everywhere and 
i?3 = 00. Wc now allow a (5-pulse to be incident from the 
left lead and calculate the transmission amplitude for the 
pulse to exit at the right lead. (In the bosonic language, 
the conductance is just proportional to the transmission 
amplitude). We know that the pulse broadens when re- 
flected from a resistive region of finite width and that the 
width of the reflected pulse is twice the width of the re- 
sistive region For simplicity, let us choose the pulse 
to be a (5-function. In order to maintain the pulse as 
a i5-function even after scattering from the resistive re- 
gion, we choose the resistivity profile for R2 to also be 
a (5- function given by r2{x2) = R2 5{x2)- We have seen 
earlier in Sec. HI] that doing this does not alter the fi- 
nal expression for the conductance as long as the total 
resistance of the patch remains the same. 

Taking the limit i?3 cx3 is like attaching stubs on ei- 
ther sides of the resistance R2 and connecting it to Fermi 



The total transmission amplitude for this system is a 
sum over transmission amplitudes for all possible paths. 
In addition to the two paths shown in Fig. [3] (i), there 
are infinitely many paths in which the pulse starting in 
wire 1 ends in wire 2. These paths are characterized by 
multiple reflections of the pulse between the junction- 
and the free end of the wire 3. One such path is shown 
in Fig.|9](ii). The path 1 ^ (3^ 3)" 2 from wire 1 
to wire 2 with n reflections between junction- T' and the 
free end of the wire 3, has a transmission amplitude of 
M-P23 (-^^■P33)" A/-P31 (this is non-zero only when O-p =/= 
0). Hence the total transmission amplitude tp{l 2) 
(when 9-p 7^ 0) turns out to be 

Ml ^2) = Mp2i + ^^^2^. (26) 

i - Mp33 

The DC conductance of such a system connected to leads 
on the two sides, is given by adc = (e^/ft.) t-p{l — >■ 
2). Using the fact that the matrix M-p describing the 
junction is parameterized by a single variable 0p, one 
can show that (i) adc = e^/ft, for Mp lying in the Mi 
class, and (ii) adc — for Mp lying in the M2 class, 
irrespective of the angle 6p as long as 9p 7^ 0. For the Mi 
class, 9p = means that the three wires are disconnected 



8 



and it is easy to see that (Jdc = 0, while for the M2 class, 
O-p = means that the stub is disconnected from the wire 
and the calculation gives the expected result, adc = e^/h. 

It is easy to see by this method that the effective resis- 
tance of a system having two or more scattercrs in scries 
is the sum of the individual resistances. In other words, 
if ti, t2, ... are the transmission probabilities of the dif- 
ferent scatterers in series, the quantities — l)'s add 
up to give 1/i — 1, where t is the transmission probabil- 
ity for the system (see pp. 63-64 of Ref. 0). Now, the 
model in Fig. [8] can be looked upon as a system with 
three scattering centers, i.e., a resistance R2 in the mid- 
dle and two y-j unctions on the two sides of i?2- When 
either of the two junctions {V) is described by an M2- 
matrix (with 9^0), we get = cx3 since the resistance 
of the scatterer at V is infinite. When both the junctions 
are described by A/i-matrix (with the constraints that 
9c ^0 and 617^ ^ 0), we get Rn = i?2. 



VII. SPIN-1/2 ELECTRONS 

Let us briefly discuss an extension of our results to the 
realistic case of spin-1/2 electrons. In one dimension, it 
is known that interactions between electrons lead to the 
phenomenon of spin-charge separation A bosonic 

description of the system begins with fields for spin-up 
and spin-down electrons denoted by (f>^ and (f>^ respec- 
tively. The fields for the charge and spin modes are then 
given by 



and (0t - M- 



(27) 



The system decouples in terms of these fields even when 
density-density interactions are introduced between the 
electrons. The Lagrangian is similar in form to the one 
in Eq. ([T]) except that there are two sets of parameters 
denoted by {Kc, v^) and {Kg, Vg) for the charge and spin 
fields respectively. Namely, 



/oo 
dx 
-00 



ZVcKc 2,Kc 
+ 77^{dt<l^s? - ^{d.<i^s)\ (28) 

where we have ignored a term involving the cosine of 
the field (j}s arising from fourth-order fermionic terms like 
(^jV'j,)^ [3- For an inhomogeneous system, the param- 
eters Ka and Va are generally functions of x. However, 
for a system with SU{2) rotational invariance. Kg = 1. 
The charge and spin currents are given by 



7 c = ^V2dt 



3c = —j= V ^ Utipc, 

a/tt 



respectively. We note that jc is invariant under SU{2) 
rotations, while js is only invariant under U{1) rotations 
about the z-axis. The simplest way of introducing re- 
sistance in this theory would be to postulate a Rayleigh 
dissipation function for the charge current given by 



J" = 



dx r 



e 

2^ 



dx 2r {dtcj)cY. (30) 



We can now use this dissipative function to calculate the 
charge conductance in a variety of systems as in the spin- 
less case. In principle, one can also introduce a dissipa- 
tive function for the spin current; however such a function 
would have both a term quadratic in the field 0^ as well 
as a cosine of (j)s [13 • This makes it difficult to analyze 
the corresponding equations of motion. 

In analogy with Eqs. (jlOmi|) . we can compute the two- 
point correlation function of the charge current jc to 
find the charge conductance Gc for a three-wire junc- 
tion. The presence of the factors of in Eq. (|29|) and 2 
in Eq. (j30[) as compared to the corresponding expressions 
for the spinless case implies that Gc will be given by an 
expression similar to Eq. (flS)) . except for some factors of 
2. Namely, we will have 



Gc 



2e^K, 



[1 + AI + Kcw{t - M){1 + — R)]-i 

TT TT 



X [1-M], 



(31) 



and js 



V2 9t 



(29) 



where the current splitting matrix M and the resistance 
matrix IR are defined as before, and Kcw denotes the Lut- 
tinger parameter for the charge field in the wire regions 
given by < Xi < Ln. 

One can actually define two conductances, Gc and 
Gs, which govern the amounts of charge and spin cur- 
rents which flow when the corresponding voltage bi- 
ases or chemical potential differences are applied be- 
tween different leads, namely, eAVc = A/i^ = A/i^ and 
e/SVs ~ A/i-|- = — A/xj, for driving charge and spin cur- 
rents respectively. In the absence of resistances, it has 
been shown in Ref. [l^j that for a junction of three quan- 
tum wires, the RG ffows resulting from the interactions 
between the electrons generally take the system to a fixed 
point at which the charge and spin conductances, Goc 
and Gos, are different from each other. [We note here 
that the RG fiows occur entirely within the interacting 
regions < Xi < Ln in Fig. [51 and that no renormal- 
ization occurs in the dissipative regions since there is no 
interaction in those regions]. 

Let us now introduce resistances in the system. If we 
then use a relation like Eq. ([TT]) to find Ga from Goa 
and R, where a = c, s, we expect to find that Gc and Gs 
would also not be equal to each other since Gqc and Gqs 
differ from each other. Note that this is purely an effect 
of interactions between the electrons; for non-interacting 
electrons, we would have Gqc = Gos and therefore Gc = 
Gs. 
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VIII. DISCUSSION 



We have developed a formalism which allows us to 
study the effect of resistances in a quantum wire using the 
technique of bosonization. The analysis can be extended 
to a system of three wires by using a current splitting 
matrix M to describe the junction. It is known that 
there are two classes of such matrices which are called 
Ml and M2. We have calculated the conductance of a 
three-wire system using both Green's function and scat- 
tering approaches. We have then examined the power 
dissipated by the system as a function of the matrix Ad 
and the voltages applied in the three leads. For the Mi 
class, both the conductance and the dissipated power de- 
pend on the value of the interaction parameter K near 
the junction, while in the M2 class, the conductance and 
power arc independent of K . Next, by putting together 
two junctions of three wires, we have studied the effective 
conductance of a system of two wires in parallel. This is 
found to depend in a highly non-trivial way on the M 
matrices at the two junctions, the parameter if, and the 
resistances in the two wires. In some cases, the effective 
resistance is infinite, while in other cases, it is finite but 
depends on K. In only one special case do we obtain the 
classical result for two parallel resistances. For the case 
in which one of the two resistances is infinitely large, we 
have provided an intuitive way of calculating the effec- 



tive conductance by summing over all the paths that an 
electron can take in going from one lead to the other. 
This method also shows that three-wire junctions with 
matrices Mi and M2 behave quite differently from each 
other when the resistance in one of the wires is taken 
to be infinitely large. Finally, we have generalized our 
results to the case of interacting spin-1/2 electrons. We 
have argued that the charge and spin conductances will 
generally be different from each other due to RG flows 
induced by interactions between the electrons. 

The formalism discussed in this paper is well suited 
for studying systems with interacting electrons in which 
the resistances are phase incoherent. It would be useful 
to develop a more general method which can deal with 
partially coherent resistances. It would also be useful to 
have a more complete treatment of systems with spin- 
1/2 electrons in which the charge and spin resistivities 
are position dependent and different from each other, in 
the spirit of Ref . [Hi . 
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